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We perform a many-body investigation of the spin dephasing in n-typed InAs quantum wells 
under moderate magnetic fields in the Voigt configuration by constructing and numerically solving 
the kinetic Bloch equations. We obtain the spin dephasing time due to the Rashba effect together 
with the spin conserving scattering such as the electron-phonon, the electron-nonmagnetic impu- 
rity as well as the electron-electron Coulomb scattering. By varying the initial spin polarization, 
temperature, impurity density, applied magnetic field and the interface electric field, we are able to 
study the spin dephasing time under various conditions. For the electron density and quantum well 
width we study, the many-body effect dominates the spin dephasing. Moreover, we find an anoma- 
lous resonance peak in the spin dephasing time for high initial spin polarization under moderate 
magnetic fields. 

PACS numbers: 71.10.-w, 67.57.Lm, 72.25.Rb, 73.61.Ey 



I. INTRODUCTION 

Almost all of the current semiconductor devices are 
based on manipulating electronic charges. The aris- 
ing field of spintronics proposes to use the spin de- 
gree of freedom of electrons in place of/in addition 
to the charge degree of freedom for device applica- 
tions in order to add new features and functionalities 
to semiconductors devices. 1,2,3 The hope of the real- 
ization of the proposed spintronic devices is supported 
by the resent development of ultrafast nonlinear optical 
experiments 4,5,6,7,8,9,10,11,12,13,14,15,16,17 where long spin 
dephasing time (> 100 ns) is reported. 

The functionalities of the semiconductor spintronic 
devices rely on the manipulation of the spin co- 
herence. In order to realize these devices, one 
needs to thoroughly understand the spin dephas- 
ing mechanisms which tend to destroy the spin co- 
herence. Historically, three spin dephasing mecha- 
nisms have been proposed in semiconductors: 18,19 the 
Elliot- Yafet (EY) mechanism, 20,21 the D'yakonov-Perel' 
(DP) mechanism, 22 and the Bir-Aronov-Pikus (BAP) 
mechanism. 23 All of the three mechanisms are either due 
to the spin flip (SF) scattering or are treated as effective 
SF scattering. The spin dephasing times of these mech- 
anisms for low polarized system are calculated in the 
framework of single particle approximation. 18 In addi- 
tional to these single particle spin dephasing mechanisms, 
three years ago Wu proposed a many-body spin dephas- 
ing mechanism which has long been overlooked in the 
literature. This mechanism is caused by irreversibly dis- 
rupting the phases between spin dipoles due to the inho- 
mogeneous broadening together with the spin-conserving 
(SC) scattering, 24,25 ' 26,27 ' 28 and is therefore a many-body 
effect. The inhomogeneous broadening can be introduced 
by the energy dependence of g-factor 24,25,29 and/or the 
momentum k-dependence of the DP term. 25,26,2 Our re- 
cent works further show that this mechanism also plays 
an important role in the spin dephasing during the spin 



transport. 30,31 

Very recently we performed a systematic 
investigation 32,33 of the spin dephasing due to the 
DP effect in n-typed GaAs (100) quantum wells for high 
temperatures (> 120 K) under magnetic fields in the 
Voigt configuration by constructing and numerically 
solving the kinetic Bloch equations. 24,25,26,27,28,34 In 
these studies, we include all the SC scattering such as 
the electron-phonon, the electron-nonmagnetic impurity 
as well as the electron-electron Coulomb scattering and 
investigate the spin dephasing under various conditions. 
The dephasing obtained from our theory contains both 
the single-particle dephasing caused by the effective SF 
scattering first proposed by D'yakonov and Perel', 22 
and the many-body one due to the inhomogeneous 
broadening provided by the DP term. We show that 
for the electron densities we studied, the spin dephasing 
rate is dominated by the many-body effect. Moreover, 
as we include the electron-electron Coulomb scattering, 
we are able to investigate the spin dephasing with extra 
large spin polarization (up to 100 %) which has not 
been discussed both theoretically and experimentally. 
We find that under moderate magnetic fields, the SDT 
increases dramatically with the initial spin polarization. 
For example, the SDT of a impurity free sample gets 
an increase of more than one order of magnitude when 
the initial spin polarization rises from about to 
about 100 % at low temperature. 33 The initial-spin- 
polarization dependence of the spin dephasing becomes 
more interesting when the magnetic field is increased 
to a few tens tesla where the SDT no longer increases 
monotonically with the initial spin polarization but 
shows an anomalous resonance peak versus the initial 
spin polarization. 32 The dramatic increase and the 
anomalous resonance of SDT in the high spin polariza- 
tion region is found to be due to the first order of the 
electron-electron interaction, i.e., the Hartree-Fock (HF) 
contribution which provides an effective magnetic field 
that can reduce the spin dephasing and result in a fast 
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increase of the SDT. Moreover, under right condition, 
the HF term, the applied magnetic field as well as the 
DP term can reach to a resonance, and thus forms the 
anomalous peak. Due to the small Lande g-factor in 
GaAs, the resonance condition can only be achieved 
under very high magnetic fields. 

In this paper, we apply the kinetic theory to study the 
spin dephasing in the n-typed InAs QW for high tem- 
peratures where the DP term is the leading dephasing 
mechanism. In QW system, the DP term is composed 
of the Drcsselhaus term 35 and the Rashba term. 36 ' 37 The 
Dresselhaus term is due to the lack of inversion symmetry 
in the zinc-blende crystal Brillouin zone and is sometimes 
referred to as bulk inversion asymmetry (BIA) term. 
Whereas the Rashba term appears if the self-consistent 
potential within a QW is asymmetric along the growth 
direction and is therefore referred to as structure inver- 
sion asymmetry (SIA) contribution. For QW's of wide 
band-gap semiconductors such as GaAs, the Dresselhaus 
term is the main spin dephasing mechanism. Whereas for 
QW's composed of narrow band-gap semiconductors such 
as InAs in the present case, the Rashba term is dominant. 
As the Rashba term is proportional to the interface elec- 
tric field of the QW, therefore, the spin dephasing in 
the InAs QW can be manipulated through applying an 
electric field perpendicular to the QW. Moreover, as the 
Lande g-factor in InAs is very large (g = 15 compared to 
0.44 of GaAs), one expects to achieve the resonance con- 
dition under a moderate magnetic field. We organize the 
paper as follows: We present our model and the kinetic 
equations in Sec. II. Then in In Sec. III(A) we investigate 
how the SDT changes with the variation of the initial 
spin polarization. The temperature dependence of the 
SDT under different spin polarization is discussed in de- 
tail in Sec. III(B), where we also highlight the difference 
between the present many-body theory and the earlier 
simplified theory. In Sec. III(C) we show the magnetic 
field dependence of the SDT. Finally we discuss how the 
interface electric field affects the SDT. We present the 
conclusion and summary in Sec. IV. 



II. KINETIC EQUATIONS 

We start our investigation from an ra-doped (100) InAs 
QW with well width a. The growth direction is assumed 
to be z-axis. A moderate magnetic field B is applied 



along the x axis. Due to the confinement of the QW, 
the momentum states along z axis are quantized. There- 
fore the electron states are characterized by a subband 
index n and a two dimensional wave vector k = (k x , k y ) 
together with a spin index a. In the present paper, we 
choose the electron density so that only the lowest sub- 
band is populated and the transition to the upper sub- 
bands is unimportant. Therefore, one only needs to con- 
sider the lowest subband. With the DP term (specifically 
the Rashba term) included, the Hamiltonian of the elec- 
trons in the QW takes the form: 



H 



\i<7<7 



;| £k +[ fl/iB B + h(k)]-^ }ck<W+ffj. (1) 



Here £k = k 2 /2m* is the energy of electron with wavevec- 
tor k and effective mass m*. u are the Pauli matrices. 
The Rashba term h(k) can be written as 

h x (k) = ak y , h y (k) = -ak x , h z (k) = 0. (2) 

In these equations, a is proportional to the interface elec- 
tric field E z along the growth direction: 



a = a eE z , 



(3) 



with the coefficient a being inversely proportional to 
the energy gap and the effective mass. 38 The interac- 
tion Hamiltonian Hi is composed of Coulomb interaction 
H ee , electron-phonon interaction H p h, as well as electron- 
impurity scattering Hi. Their expressions can be found 
in textbooks. 39,40 

We construct the kinetic Bloch equations by the 
nonequilibrium Green function method 39 as follows: 



Pk,aa' — Pk : erer'lcoh + P\l.<j<j' Iscatt 



(4) 



Here p k represents the single particle density matrix. 
The diagonal elements describe the electron distribu- 
tion functions p k<J<J = ,/kcr- The off-diagonal elements 
p k i_i = /?k describe the inter-spin-band polarizations 

(coherence) of the spin coherence. 34 Note that p k 



11 = 

2 2 



= p k . Therefore, / k ±i an d Pk are the quantities 
to be determined from Bloch equations. 

The coherent part of the equation of motion for the 
electron distribution function and the spin coherence are 
given by 



J 



dt 



co h = -2a{[gn B B + /i x (k)]Imp k + h y (k)Rcp k } + 4CTim^V q p k+q/ o k , 



(5) 



d 1 

-^f Icoh = 2-feMB£ + «M k )+M k )K/ki - /k-j)+*$Zlq[(/k+q* ~ /k+q-1 Vk ~ />k+q(/ k I ~fk-$)], ( 6 ) 
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respectively, where Vq = Aire 2 /[no(q + qo)] is the 2D 
Coulomb matrix element under static screening, qo = 
(e 2 m*/Ko) J2a /k=o,cr and Ko is the static dielectric con- 
stant. The first term on the right hand side (RHS) of 
Eq. (5) describes the spin precession of electrons under 
the magnetic field B as well as the effective magnetic field 
h(k) due to the Rashba effect. The scattering terms of 
the electron distribution function and the spin coherence 
are given by Eqs. (Al) and (A2) in the Appendix. 
The initial conditions are taken at t = as: 

Pk\t=o = (7) 

/ka|t=o = l/{exp[(e k - n a )/k B T\ + 1} (8) 

where \i a is the chemical potential for spin a. The condi- 
tion fii gives rise to the imbalance of the electron 
densities of the two spin bands. 



III. NUMERICAL RESULTS 

The kinetic Bloch equations form a set of nonlinear 
equations. All the unknowns to be solved appear in the 
scattering terms. Specifically, the electron distribution 
function is no longer a Fermi distribution because of the 
existence of the anisotropic Rashba term h(k). This term 
in the coherent part drives the electron distribution away 
from an isotropic Fermi distribution. The scattering term 
attempts to randomize electrons in k-space. Obviously, 
both the coherent part and the scattering terms have 
to be solved self-consistently to obtain the distribution 
function and the spin coherence. 

We numerically solve the kinetic Bloch equations in 
such a self-consistent fashion to study the spin preces- 
sion between the spin-up and -down bands. We include 
electron-phonon scattering and the electron-electron in- 
teraction throughout our computation. As we concen- 
trate on the relatively high temperature regime in the 
present study, for electron-phonon scattering we only 
need to include electron-LO phonon scattering. Electron- 
impurity scattering is sometimes excluded. As discussed 
in the previous paper, 34 ' 41 irreversible spin dephasing can 
be well defined by the slope of the envelope of the inco- 
herently summed spin coherence p(t) = |pk| • The 
material parameters of InAs for our calculation are tab- 
ulated in Table I. 42 The method of the numerical calcu- 
lation has been laid out in detail in the previous paper 
on the DP mechanism in 3D systems. 24 The difference 
is that here we are able to get the results quantitatively 
in stead of only qualitatively as in the previous 3D case, 
thanks to the smaller dimension in the momentum space. 
Our main results are plotted in Figs. 1 to 6. In these cal- 
culations the total electron density N e and the applied 
magnetic field B are chose to be 4 x 10 10 cm -2 and 1.5 T 
respectively unless otherwise specified. 
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TABLE I: Parameters used in the numerical calculations 



A. Spin polarization dependence of the spin 
dephasing time 

We first study the spin polarization dependence of the 
SDT. As our theory is a many-body theory and we in- 
clude all the scattering, especially the Coulomb scatter- 
ing, in our calculation, we are able to calculate the SDT 
with large spin polarization. 

In Fig. 1, SDT r is plotted against the initial spin 
polarization P with A, = [Fig. 1(a)] and A, = 0.1 N e 
[Fig. 1(b)] at different temperatures. The most striking 
feature of the impurity-free case is the huge anomalous 
peaks of the SDT in low temperatures. For T = 120 
K, the peak value of the SDT is about 6 times higher 
than that of low initial spin polarization. It is also seen 
from the figure that the anomalous peak is reduced with 
the increase of temperature and the peak shifts to higher 
polarization. For T > 200 K there is no anomalous peak. 

The anomalous peak in the t-P curve in low tem- 
perature region originates from the electron-electron in- 
teraction, specifically the Hartree-Fock (HF) self-energy 
[i.e., the last terms in the Eq. (5) and (6)]. If one 
removes the HF term, the anomalous peak as well as 
the large increase of SDT disappears. It is pointed out 
in our previous paper that although the HF term itself 
does not contribute to the spin dephasing directly, 25 ' 28 
it can alert the motion of the electrons as it behaves as 
an effective mag netic field B HF (k). Therefore, the HF 
term can affect the spin dephasing by combining with the 
DP term. For small spin polarization as commonly dis- 
cussed in the literature, the contribution of the HF term 
is marginal. However, when the polarization gets higher, 
the HF contribution becomes larger. Especially the effec- 
tive magnetic field formed by the HF term contains a lon- 
gitudinal component [-Bp F (k)] which can effectively re- 
duce the "detuning" of the spin-up and -down electrons, 
and thus strongly reduces the spin dephasing, therefore 
the SDT increases with initial spin polarization. 33 More- 
over, besides the initial polarization, p k and therefore 
£?H^(k) arc a } so a ff ec t ec j by the applied magnetic field. 
With higher magnetic field, both gets larger. Under the 
high magnetic field and when the initial spin polariza- 
tion reaches to a right value, the effective magnetic field 
B HF (k) may reach the magnitude comparable to the 
contribution from the DP term as well as the applied 
magnetic field in the coherent parts of the Bloch equa- 
tions and reduces the anisotropic caused by the DP term. 
Therefore, one gets much longer SDT. However, if one 
further increases the initial polarization, the HF term 
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the temperature is high enough, even largest polariza- 
tion P — 100 % cannot make the HF term to reach the 
resonance condition. Therefore, the peak disappears. 

The t-P curve is much different when the impurities 
are introduced. It is seen from Fig. 1(b) that, when the 
density of impurity is large, say Ni = 0.1N e , the fast 
rise in t-P curve remains. Nevertheless the increase is 
smaller than the corresponding one when the impurities 
are absent. In addition to the reduction of the rise in t-P 
curves, the impurities destroy the anomaly too. One can 
easily see that, with the impurity level Ni = 0.1N e , when 
the temperature is 120 K, the anomalous peak is much 
flatter than the impurity free sample, while for all other 
temperatures we study, the SDT increases uniquely with 
the polarization. 
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FIG. 1: Spin dephasing time r versus the initial spin polar- 
ization P with different impurity concentration and different 
temperatures. The impurity densities in (a) and (b) are and 
0.1 A^e respectively. The lines are plotted for the aid of eyes. 



exceeds the resonance condition. As the result, the SDT 
decreases. Therefore, one gets the anomalous peak which 
is similar to the resonance effect. It is noted at, as both 
the DP term and the HF term are fc-dependent, the res- 
onance is broadened. 

For high temperatures the HF term is smaller. In order 
to reach the resonance, one needs to go to higher polar- 
ization. Therefore, as shown in the figure the anomalous 
peak shifts to the higher polarization. However, when 



FIG. 2: Spin dephasing time r versus the initial spin po- 
larization P for a InAs QW with different impurity levels. 
Circle (•): Ni = 0; Diamond (♦): AT* = 0.05N e ; Square (■): 
N = O.lAe; Up triangle (A): Ni = 0.2AT e . The lines are 
plotted for the aid of eyes. 

To further reveal the contribution of the impurity to 
the dephasing under different conditions, we plot the 
SDT as a function of the polarization for different im- 
purity levels at T = 120 K in Fig. 2. The figure clearly 
shows that the impurity tends to remove the anomalous 
peak and to shift the peak to the larger initial spin po- 
larization. This is because that the impurity reduces the 
HF term and therefore the resonance effect is also re- 
duced. Hence, in order to reach the maximum resonance, 
one has to increase the initial spin polarization. Conse- 
quently, the peak shifts to larger P. Whereas when Ni 
is raised to 0.2N e , the HF term is reduced too much to 
form a peak. 
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B. The temperature dependence of the spin 
dephasing time 

Above we discussed the dependence of spin dephas- 
ing on initial spin polarization for different temperatures. 
Now we turn to the temperature dependence of the SDT 
under different initial spin polarizations. ^From Fig. 1(a) 
and (b) in Sec. III(B), one can see that for small polariza- 
tion, the SDT increases with the temperature. Whereas 
in high polarized region, the SDT decreases with the 
temperature. For moderate polarization, the tempera- 
ture dependence is too complicated to be described by a 
monotonic function of temperature. 

To see more detail of how the spin dephasing depends 
on the temperature, we replot the SDT shown in Fig. 1 
as a function of the temperature for different impurity 
levels and different spin polarizations in Fig. 3(a) and 
(b). It is seen from the figure that similar to the case of 
GaAs where the DP term is composed of the Dresselhaus 
term, 33 for Rashba term and for low spin polarization the 
SDT also increases with the temperature for all impurity 
levels. This property is again opposite to the results of 
earlier simplified treatments of the DP effect, where it 
was predicted that the spin lifetime decreases with the 
increase of temperature in the 2D system. 43,44 The SDT 
based on the simplified model is given by 27 ' 43 ' 45 



in which 



and 



1 _S^°dE k (f H -f k _ i )r(k) 
r f °° dE k (f k , - f k _ h ) ■ 



T{k) = 2T 1 (k)(a Ek) 2 



/■2-7T 

r- 1 (k)= a(E k ,9)[l-cos(nd)]d9 . (11) 
Jo 

<j(E kl 9) stands for scattering cross-section. For compari- 
son, we plot the SDT predicated by the earlier model and 
by our present many-body theory in the inset of Fig. 3(a). 
^From the inset one can see that the SDT predicated by 
the earlier model is about one order of magnitude larger 
than the one predicated by our theory. In the mean time, 
the SDT of the earlier mode drops dramatically with the 
increase of the temperature. Nevertheless, in our many- 
body treatment, it rises slightly with the temperature. 

The giant difference between two models lies on the 
fact that the earlier simplified model is based on the 
single particle picture which does not count for the de- 
phasing due to the inhomogeneous broadening inher- 
ited in the Rashba term, which is exactly the result of 
many body effect. 24 ' 25, 26,27,28 By comparing the theoret- 
ical SDT predicated by the two models, we can see that 
the spin dephasing due to the inhomogeneous broaden- 
ing is much more important. In the case we calculated, 
the spin dephasing is dominated by the inhomogeneous 
broadening. Therefore, it is easy to understand why the 



(9) 



(10) 



earlier simplified treatment of the DP mechanism gives 
much slower spin dephasing. Although there is no exper- 
iment results for InAs, we have shown that for GaAs the 
prediction of our many-body theory agrees both quanti- 
tatively and qualitatively with the experiment results. 33 
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FIG. 3: Spin dephasing time r versus the temperature T for 
InAs QW's with spin polarization P = 2.5 % (a) and P — 
75 % (b) under two different impurity levels. Circle (•): Ni = 
0; Diamond (♦): Ni = 0.1iV e . The lines are plotted for the 
aid of eyes. The SDT predicated by the simplified treatment 
of DP term (solid curve) and our model (circle) for Ni = is 
plotted in the inset of (a) for comparison. 

The temperature dependence of the SDT can be easily 
understood when the spin dephasing due to inhomoge- 
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neous broadening is taken into account: When the tem- 
perature increases, the inhomogcncous broadening is re- 
duced as the electrons are distributed to the wider Re- 
states. As a result, the number of electron occupation 
on each k state is reduced. It is further noted that this 
reduction is mild as function of the temperature. There- 
fore, the temperature dependence is quit mild unless it 
is within the regime of anomalous peak. 

In the region where HF term is important, in addition 
to the above mentioned two effects of the temperature 
acting on the spin dephasing, the temperature depen- 
dence of HF term should also be taken into account when 
we study the temperature dependence of the spin dephas- 
ing. In high spin polarization region, the SDT decreases 
with the temperature. However, in the moderately polar- 
ized region, the temperature dependence of SDT due to 
the combination of these three effects is too complicate 
to be described by a monotonic function. We replot the 
SDT as a function of the temperature in Fig. 3(b), for the 
high polarization P — 75%, which is near the anomalous 
peak shown in Fig. 1(a). We can see that due to the re- 
duction of the HF term, the resonance is removed and the 
SDT drops dramatically with the increase of the temper- 
ature in the impurities free sample. While for the system 
with impurity concentration Ni = 0.1N e , the HF term is 
less important than that in the impurities free sample, 
and the SDT is less sensitive to the temperature. 
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FIG. 4: Spin dephasing time r versus the applied magnetic 
field for InAs QW's for different spin polarizations and differ- 
ent impurity levels. Solid curve with dots: Ni = 0, P = 2.5 %; 
Solid curve with diamonds: N = 0.1N e , P = 2.5 %; Dashed 
curve with dots: N — 0, P = 75 %; Dashed curve with 
diamonds: N =0.1N e ,P = 75 %. 



C. Magnetic field dependence of the spin dephasing 

We now investigate the magnetic field dependence of 
the spin dephasing. In Fig. 4, we plot the SDT versus 
the applied magnetic field for different impurity levels 
and different spin polarizations. It is seen that for all 
the cases we study, the SDT increases with the magnetic 
field. This is because in the presence of a magnetic field, 
the electron spins undergo a Larmor precession around 
the magnetic field. This precession suppresses the preces- 
sion about the effective magnetic field h(k). 18,29 There- 
fore the SDT increases with the magnetic field. It is 
pointed out that in 3D electron gas, the magnetic field 
also forces electrons to precess around it. This preces- 
sion introduces additional symmetry in the momentum 
space that limits the k-space available to the DP effect 
which is anisotropic in it. 18 ' 24 ' 29 This can further reduce 
the spin dephasing. However, it is expected that this ef- 
fect in the 2D case is less effective than the 3D case as 
in z-direction the momentum is quantized and the mo- 
mentum precession around the magnetic field should be 
suppressed. 

In additional to the above mentioned effect of the mag- 
netic field on spin dephasing, one can further see from 
Fig. 4, that for large polarization, the magnetic field also 
enhances the HF term. As we mentioned before, for large 
polarization, the contribution from the HF term is im- 
portant. Increase of the HF term serves as additional 
magnetic field which further suppresses the effect of the 
DP term h(k), and therefore results in a faster rise in 
the t-B curve in the region of B < 4 T. When the ap- 
plied magnetic field exceeds 4 T, the increase of HF term 
saturates, thus the slop of the t-B curves in the region 
of B > 4 T is reduced to that in the low polarization 
region. To reveal more concrete about the combining 
effect of the magnetic field and the HF term on spin de- 
phasing, we plot the SDT as a function of polarization 
in Fig. 5. It is shown that the rise in the t-P curve 
increases with the magnetic field. Moreover, the posi- 
tion of the peak in t-P shifts to a larger polarization. 
This is understood that, it needs a larger HF term, and 
hence a larger spin polarization, to achieve the resonance 
condition when the magnetic field increases. When the 
magnetic field is raised to 4 T, it is no longer possible 
to form the resonance for all of the polarization. As a 
result the SDT increases uniquely with the polarization 
and there is no peak in the t-P curve. 



D. Interface electric field dependence of spin 
dephasing 

We now investigate how the interface electric field af- 
fects the spin dephasing. In Fig. 6 we plot the SDT as 
a function of the initial spin polarization for three differ- 
ent interface electric fields. It can be seen from the fig- 
ure that when the interface electric field decreases from 
2 x 10 4 V/cm to 5 x 10 3 V/cm, the SDT increases about 
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FIG. 5: Spin dephasing time r versus the polarization for 
InAs QW's at different magnetic field. Circle (•): B=1.5 T; 
Diamond (♦): B=2 T; Square (■): B=4 T. The lines are 
plotted for the aid of eyes. 



16 folds. It is understood that when the interface elec- 
tric field increases, the Rashba effect is enhanced, and 
consequently the spin dephasing is also enhanced. More- 
over, for fixed initial spin polarization, the HF term is 
reduced when E z is increased. As a result, the interface 
electric field E z also changes the anomalous peak in the 
t-P curve as in order to achieve the resonance condition, 
one has to go to higher polarization in order to get large 
enough HF term. Consequently the resonance peak is 
smoothed and the position moves to higher initial spin 
polarization region when E z increases, which is shown in 
the Fig. 6. 



IV. CONCLUSION 

In conclusion, we have performed a systematic in- 
vestigation of the Rashba effect on the spin dephasing 
of n-typed InAs QW's under moderate magnetic fields 
in Voigt configuration. Based on the noncquilibrium 
Green's function theory, we derived a set of kinetic Bloch 
equations for a two-spin-band model. This model in- 
cludes the electron-phonon, electron-impurity scattering 
as well as the electron-electron interaction. By numer- 
ically solving the kinetic Bloch equations, we study the 
time evolution of electron densities in each spin band 
and the spin coherence - the correlation between spin-up 
and -down bands. The spin dephasing time is calculated 
from the slope of the envelope of the time evolution of 
the incoherently summed spin coherence. We therefore 



FIG. 6: Spin dephasing time r versus the initial spin polar- 
ization at T = 120 K, B = 1.5 T, Ni = for three different 
interface electric fields. 

are able to study in detail how this dephasing time is af- 
fected by spin polarization, temperature, impurity level, 
magnetic field and interface electric field. Differing from 
the earlier studies on spin dephasing based on the single 
particle model which only considers the effective SF scat- 
tering, our theory also takes account of the contribution 
of many-body effect on the spin dephasing. 24 ' 25 ' 26 ' 27 ' 28 In 
fact, for the n-typed semiconductors and the spin polar- 
ization studied in the experiments, this many-body de- 
phasing effect is even more important than the effective 
SF scattering as it is one order of magnitude larger than 
the later. Equally remarkable is that, as we include all 
the scattering, especially the Coulomb scattering in our 
many-body theory, now we are able to calculate the spin 
dephasing with extra large (up to 100 %) initial spin po- 
larization which is unable to be calculated by the former 
single particle theory. 

It is discovered that the SDT increases with the ini- 
tial spin polarization. Moreover, for low impurity level 
and low temperature, there is a giant anomalous resonant 
peak in the curve of the SDT versus the initial spin polar- 
ization. This resonant peak moves to high spin polariza- 
tion and its magnitude is fast reduced (enhanced) until 
the whole resonance disappears if one increases the im- 
purity density, the temperature and/or the interface elec- 
tric field (the magnetic field). It is discovered that this 
anomalous resonance peak originates from the HF contri- 
bution of the electron-electron Coulomb interaction. Un- 
der the right spin polarization, the contribution of HF 
term may reach the magnitude comparable to the con- 
tribution of the DP term as well as the magnetic field in 
the coherent part of the Bloch equation and reduces the 
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anisotropy caused by the Rashba effect — consequently re- 
duces the spin dephasing. As the resonance is the com- 
bined effects of the HF term, the Rashba term and the 
magnetic field, the magnitude and position of the reso- 
nance peak are affected by all the factors that can alert 
the magnitude of the HF term, such as temperature, im- 
purity scattering, magnetic field as well as the interface 
electric field: For a given impurity concentration, when 
the temperature increases, the HF term reduces. Conse- 
quently the t-P curve is smoothed and the peak position 
is moved to higher spin polarization; For impurity free 
sample, if the temperature is raised to 200 K, the HF 
term is reduced too much to form a resonance and the 
anomalous peak disappears; The same situation happens 
when the impurity level increases at a given temperature 
as the scattering also lowers the HF term. When the im- 
purity level is raised to 0.2iV e there is no resonance in 
the temperature region we studied; While the increase 
of the magnetic field enhances the HF term and results 
in a faster increase of the SDT as well as a higher reso- 
nant peak in t-P curve. Moreover, as the magnetic field 
becomes larger, it needs a larger HF term and hence a 
larger polarization in order to achieve the resonant condi- 
tion. Therefore the peak position is also moved to higher 
polarization; When the interface electric field increases, 
the HF term is reduced. Therefore the resonance peak 
in t-P becomes flatter and its position moves to higher 
spin polarization. 

For low spin polarized regime, the SDT increases when 
the temperature rises. This is contrary to the result 
of earlier simplified single-particle calculation where the 
SDT always decreases with the increase of the temper- 
ature. Moreover, the SDT predicted by our many-body 
calculation is one order of magnitude faster than the ear- 
lier result. The physics of this feature is due to the addi- 
tional many-body spin dephasing channel due to the in- 
homogeneous broadening provided by the Rashba term, 
which by combining with the SC scattering also causes 
spin dephasing. In the situation we studied, the spin de- 
phasing is dominated by the many-body dephasing effect. 
With the increase of the temperature, the inhomogeneous 
broadening reduces and the SDT increases. 

In high spin polarization region, the contribution of 



the HF term should be taken into consideration. This 
brings more complication in the study of the spin de- 
phasing. Usually the SDT can not be described by a 
monotonic function of the temperature and the impu- 
rity concentration in high polarization regime. For po- 
larization near the resonance peak in low temperature 
and impurities free samples, the SDT decreases dramat- 
ically with the temperature as the resonance is removed 
when the temperature increases. Whereas when the im- 
purity concentration is 0.17V e , the SDT is less sensitive 
to the temperature. 

As the magnetic field causes the electron spins to pre- 
cess about it, this precession will suppress the precession 
about the effective magnetic field h(k) originated from 
the Rashba effect. As a result the spin dephasing is re- 
duced. The magnetic field also enhances the HF term, 
which servers as an additional magnetic field and fur- 
ther suppresses the Rashba effect. Therefore, the t-B 
curve gets a faster increase in the high polarization re- 
gion. Our calculation also shows that when the interface 
electric field increases, the SDT decreases. This is be- 
cause with the increasing of the interface electric field, 
the Rashba term is strengthened. 

In summary we have performed a thorough investiga- 
tion of the spin dephasing in n-typed InAs QWs. Many 
new features which have not been investigated both theo- 
retically and experimentally are predicted in a wide range 
of parameters. 
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APPENDIX A: APPENDIXES 



J 



dt 



The scattering terms of electron distribution functions in the Markovian limit are given by 

<9/k,<T | 

scatt - ) 

<iq z \ 

-2wNi ^2 UqS(e k - £ k _ q ) [/k<r(l - /k-q<r) - Re(p k/ 9 k _ q )] - 27T V q S ( £ k-q - £k + £k' ~ £k' 

q qk'tr' 



= \ ~ 27F ^QA^k - £k-q - ftqg^A) [Nq gx x{fka ~ /k-q CT ) + /ka(l - /k-qa) - Rc(/3 k p k _ q )] 



(1 - /k-q CT )/kcr(l - /k'crO/k'-qo-' + T^kPk-qt/k'a' - /k'-qcr') + ^k'Pk'-ql/k-qa _ /kf) 

— {k <-> k — q, k' <-» k' — q}, 



(Al) 



in which {k <-> k— q, k' <-> k' — q} stands for the same terms as in the previous {} but with the interchange k <-> k— q 
and k' <-» k' — q. The first term inside the braces on the RHS of Eq. (Al) comes from the electron-phonon interaction. 
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A stands for the different phonon modes, i.e., one longitude optical (LO) phonon mode, one longitudinal acoustic (AC) 
phonon mode due to the deformation potential, and two AC modes due to the transverse piezoelectric field. g m ^\ 

arc the matrix elements of clcctron-phonon coupling for mode A. For LO phonons, = {4naQ^Q/[y/2jl(q 2 + 

ql)]}\I(iq z )\ 2 with a = e 2 y/ fi / (2CIj j q) (k^ 1 — k^ 1 ). is the optical dielectric constant and ^lq is the LO phonon 
frequency. The form factor \I{iq z )\ 2 = tt 2 sin 2 y/[y 2 (y 2 — tt 2 ) 2 ] with y = q z a/2. N mz \ = l/[exp(f2 qgz A/&BT) — 1] is 
the Bose distribution of phonon mode A at temperature T. The second term inside the braces on the RHS of Eq. (Al) 
results from the electron-impurity scattering under the random phase approximation with Ni denoting the impurity 

concentration. £/ 2 = ^2 {AnZie 2 /[ko(q 2 + <7 2 )]} |-^(*<7z)| 2 is the electron-impurity interaction matrix element with Zi 
stands for the charge number of the impurity. Zi is assumed to be f throughout our calculation. The third term is 
the contribution of the Coulomb interaction. Similarly, the scattering parts of the spin coherence are given by 

f 

-JjT IsCatt = j 71 " ffq-J^Ok - £k-q - ^q,J,A)[pk-q(/ k l + /k- 1 ) + (/k-q| + /k-q-| - 2)p k - 2Nq qzX (pk - Pk-q)] 

+7TN t ^q<5(£k - £k- q ) [(/ k l + / k -|)Pk-q - (2 - / k _ q l - / k _ q _| )(*] 

q 

-5Z 7rV q <5 ( £k - < l " e k + ek' -£ k '- q )N/ k _ q iPk + Pk- q / k _i)(/ k 'i -/ k '-qi +/ k '-i -/ k '-q-i) 
qk' V " 

+Pk[(l - / k '|)/ k -ql + (1 - / k '-i)/ k -q-I - 2Rc(p k ,p k _ q )] - p k -q[/ k ,i (1 - / k ,_ q i) 

+(l-/k'-i)/k'- q -i-2Re(p k ,p k '- q )])J-{k^k-q,k'^k'-q} . (A2) 
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